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Present-day angle-resolved photoemission spectroscopy (ARPES) has offered a tremendous ad- 
vance in the understanding of electron energy spectra in cuprate superconductors and some related 
compounds. However, in high magnetic field, magnetic quantum oscillations at low temperatures 
indicate the existence of small electron (hole) Fermi pockets seemingly missing in ARPES of hole 
(electron) doped cuprates. Here ARPES and quantum oscillations are reconciled in the framework 
of an impurity band in the charge-transfer Mott-Hubbard insulator. 
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ARPES of cuprate superconductors [2] proved to be 
particularly instrumental in modeling the electron en- 
ergy spectrum of these charge transfer Mott-Hubbard 
insulators showing that below optimal doping the single- 
particle Fermi surface is reduced to four disconnected 
peculiar-shaped spots. While on the overdoped side of 
the phase diagram a large Fermi surface is expected, a 
few holes doped into the insulator would naturally give 
rise to four hole nodal Fermi pockets with a small area 
proportional to the doping on the underdoped side of the 
phase diagram. Another possibility is a truncation of the 
large hole Fermi surface giving rise to four nodal Fermi 
arcs due to a highly anisotropic quasiparticle life-time 
and/or a d- wave-like pseudogap. 

Remarkably, magnetic quantum oscillations (MQO) 
in kinetic and magnetic response functions of oxygen- 
ordered ortho-II YBa2Cu306.5 (YBC06.5) and of some 
other cuprates [3] revealed small electron Fermi pockets, 
rather than hole pockets or arcs seemingly in disagree- 
ment with ARPES results 0. 

A number of Fermi surface reconstructions 0] and non- 
Fermi- liquid models, including our modulated vortex lat- 
tice scenario [5| , have tempted to account for the nature 
of these unusually slow MQO. Further careful experi- 
ments have found the Zeeman splitting in MQO, which 
separates spin- up and spin-down contributions, indicat- 
ing that electrons in cuprates behave as nearly free spins, 
which rules out most of the reconstruction and non- Fermi 
liquid scenarios [fj. This observation as well as excellent 
fits of MQO data with field-independent oscillation fre- 
quencies support the view [3| that MQO are a signature 
of the true zero-field normal state. This view is also sup- 
ported by the observation of similar slow oscillations in 
electron-doped cuprates Nd2- x Ce :c Cu04 0, where due 
to their lower critical fields, the normal state is reliably 
accessed for any doping level. In contrast with the elec- 
tron pockets in hole doped cuprates, MQO reveal small 
hole pockets around certain doping level x ~ 0.165 in 
electron doped cuprates 0|. 

Until recently these two different measurement tech- 
niques, ARPES and MQO, were carried out on different 



materials. Aiming to resolve the outlined puzzle Hossain 
et al. [9] succeeded to control the surface doping and 
follow the evolution of ARPES from the overdoped to 
the underdoped regime through an in situ deposition of 
potassium atoms on cleaved YBCO. Ref. [9[ did not find 
any sign of the electron pockets in the ARPES data from 
underdoped YBCO, nor any sign of extra zone folding 
due to the kind of density wave instabilities that might 
give rise to such a Fermi surface reconstruction. Hossain 
et al. argued that if any pocket had to be postulated on 
the basis of their ARPES data, the most obvious possibil- 
ity would be that the Fermi arcs are in fact nodal Fermi 
pockets. However, these are hole, not electron pockets. 

Apparently without a detailed microscopic theory both 
ARPES and MQO data remain a mystery. Here I rec- 
oncile these two rather precise techniques by introducing 
an impurity band in the doped charged transfer Mott- 
Hubbard insulator. Small electron (hole) pockets in the 
impurity band account for MQO observed in hole (elec- 
tron) doped cuprates [3|, l7j , and for the half- moon spots 
in ARPES [3 . 

Cuprate superconductors are strongly correlated 
electron systems where the ab-initio local-density- 
approximation (LDA) fails especially in the undoped 
regime. Fortunately, adding the on-site Coulomb repul- 
sion U to the LDA analysis within the LDA+U algorithm 
[Io| or using LDA plus the tight-binding (cluster) approx- 
imation (LDA+GTB algorithm [H|), one can reproduce 
the correct magnetic ground state and the charge-transfer 
gap in parent insulators such as La2Cu04 (LCO). The 
latter and some other schemes [12j for the electronic 
structure found the charge-transfer gap also in the para- 
magnetic state of cuprates, pointing to a persistent Mott 
physics also at finite doping [l3[ . 

Different from the reconstruction models, proposed so 
far, I suggest that YBC06.5 is the Mott-Hubbard insu- 
lator where the Fermi level is pinned within the charge- 
transfer gap as in its parent insulator YBCO6.0 Q. 
This assumption is supported by experiments with ul- 
trathin insulating La2Cu04 and doped superconducting 
La2- x Sr :c Cu04 (LSCO) layers [14] and by earlier optical 
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FIG. 1: (Color online) The ortho-II phase of YBa2Cus06.5 is 
characterized by a periodic alternation of empty Cu and filled 
Cu-0 fr-axis chains doubling the unit cell in the a direction. 
Extra oxygen in the full chains gives rise to an attractive 
potential for the holes (upper panel) creating a coherent band 
inside the charge-transfer gap of the Mott-Hubbard insulator 
(lower panel). 



spectroscopies of YBCO jTH, [HI indicating that doped 
electronic states appear within the charge-transfer gap. 
These states, localized in band-tails by disorder, readily 
account for sharp quasiparticle peaks and a rapid loss of 
their intensities in some directions of the Brillouin zone 
observed in ARPES of lightly doped LSCO [17]. The 
same band-tail model explains two energy scales, their 
temperature and doping dependence, the asymmetry and 
inhomogeneity of tunneling spectra of cuprate supercon- 
ductors [18[. Pinning the Fermi level within the impu- 
rity band-tail is also compatible with the insulating-like 
low-temperature normal state resistivity and many other 
kinetic properties of lightly doped cuprates [I?], EH, Hq| • 
In the LDA picture the extra oxygen in the full chains 
of YBa2Cu306.5, Fig.JTJ splits the planar C11O2 metallic 
band in two with a gap opening at the new Brillouin zone 
boundary k x = 7r/2a, estimated from local-density ap- 
proximation calculations between 120 and 160 meV |21j . 
Quite a different band structure emerges in the charge- 
transfer insulator. Our key point is that, at variance 
with other doping levels, 0.5 (per unit cell) extra oxy- 
gen creates a coherent band within the charge-transfer 
gap rather than the localized band-tail, since YBC06.5 
is perfectly ordered, Fig.(pQ). As shown below, this coher- 
ent band accounts for the observed quantum oscillations, 
and its spectral function combined with the matrix ele- 



FIG. 2: (Color online) The impurity-band energy dispersion 
in YBa2Cu306.5 with four small electron pockets at the re- 
duced Brillouin zone boundary (dark oval spots) as described 
by Eq.® with t' = 0.7* and the Fermi level u = -2M (En- 
ergy is measured in 2t and a = 1). 



ments also explains ARPES. 

To get an insight into the coherent "impurity" band 
dispersion of YBC06.5 we employ the LDA+U band 
structure of the parent insulator YBCO6.0 [4]. In such a 
framework, the C11O2 in-plane states are found about 0.5 
eV below the Fermi level and the first electron-removal 
(i.e valence) state has a BaO — Cu chain character (see 
Fig. 5 in Ref. [4]). This valence band has its maximum 
at T point of the Brillouin zone. Spin fluctuations and 
significant electron-phonon interactions [19j could push 
the in-plane states closer to the Fermi energy, thus back 
into the first ionization energy range, due to polaronic 
level shifts missing in the LDA+U band structures. The 
edge of the upper in-plane band is also found at Y point 
in YBC06 [4(, therefore which particular orbitals form 
the valence band is not an issue in our further analysis. 
In any case we are dealing with quasi-two dimensional 
carriers. 

A single extra oxygen ion in the chain creates the 
Coulomb attractive potential for a hole, about — 2e 2 /ep, 
at a large enough distance p from the ion. Solving the 
2D Coulomb problem in the effective mass (m) approx- 
imation yields an estimate of the localized state ioniza- 
tion energy, Ei m = 8me 4 /ft 2 e 2 . Ei m is about 87 meV 
with m = 2m e measured in MQO experiments Q and 
with e = ye^ec, where e a b = 500 [22] and e c = 10 [23| 
are the in-plane and c-axis dielectric permittivities, re- 
spectively, measured in the insulating YBCO single crys- 
tals. This energy is much lower than the charge-transfer 
gap, which is about 1 eV [4j, so that the bound state is 
rather shallow. The size of the bound-state wave func- 
tion, fim(p) exp(— p/dim), is comparable to the lattice 
constant, aim = h 2 e/Ame 2 ~ 0.4 nm. 

Excess oxygen ions in YBC06.5 give rise to a poten- 
tial V(~ft), Fig.(pQ), periodic in a plane ~ft = {x,y}. It is 
"gentle" [24j in the outlined sense due to the high polar- 
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FIG. 3: (Color online) The momentum-distribution map of 
the impurity-band spectral function at the Fermi surface with 
t' 0.7t, /x -2M and 7 0.04t (here a = 1). 



izability, e ^> 1 of perovskites, so that the valence-band 
states of the parent insulator are mainly involved. Hence 
the impurity-band wave function, ipfr) can be expanded 
in a complete set of the orthogonal valence-band Wan- 
nier functions, w(r), which account for most of the cor- 
relations . These functions are atomic-like with their 
extension, ao smaller than the lattice constant, a. Thus 

V(r) = ^F(m) W (r-m), (1) 

m 

where the envelope function F(m) satisfies the following 
differential equation [25| 

[E(-i$ P ) + V(t)]F(t) = E tm F(t)- (2) 

Here E(p) is the LDU+U dispersion of the 2D valence 
band 0. 

The eigenstates of Eq.(|2|) can be expanded as 

F k (^) = ^=5>* k '"/(y-n) (3) 

since the potential V(~ft) is periodic. Here k is the quasi- 
momentum in the reduced 2D Brillouin zone, (— ir/2a < 
k x < 7r/2a, — 7r/a < k y < 7r/a), f(j} — n) are the orthog- 
onal impurity-band Wannier orbitals built of the bound- 
state wave function fim(~ft), and n are 2D position vec- 
tors of excess-oxygen ions A^ m or their projections in 
the plane, Fig.(pQ). Then the impurity-band dispersion is 
found as 

£ im (k) = -^e ik ' n t(n), (4) 

n 

where t(n) = -Jd^f(f + n)[E(-tf „) + Vtf)]ftf) 
are the impurity-band hopping integrals, which are pos- 
itive for the attractive V(~ft). It should be noted that 
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FIG. 4: (Color online) The ARPES momentum-distribution 
map of the impurity-band (upper panel), Eq.jSj) with 
(aim/a) 2 = 0.4, t' = 0.7t and 7 = 0.94t. Lower panel: the 
experimental map of the Fermi surface in YBC06.5 0] (here 
a 1). 

the impurity-band dispersion, Eq.((4]) holds for the whole 
Brillouin zone as long as we do not expand the valence- 
band energy operator E(-i^) to any finite order in 

The envelope function Eq.(j2j) works correctly to any 
power in k as long as the extension of the valence-band 
Wannier orbitals is much less than the extension of the 
impurity-band orbitals, ao <C a^ m ~ a. 

One can parametrise the impurity-band dispersion by 
keeping nearest, t and next nearest, t' , hopping integrals 
in Eq.(|4]). Thus 

E^ m (k) = -2tcos(k y a)-2t f [cos(2k x a) + cos(2k y a)], (5) 

for holes with the lower band-edge at T point, which is 
compatible with the LDA+U band structure of YBCO6.0 
[3]. Every extra oxygen donates two holes, so that the 
impurity band is almost full with holes in YBC06.5 while 
some holes are still localized by residual disorder in the 
chains and bound into preformed pairs (bipolarons) [171 - 
Il9j . The electron impurity-band dispersion, 

Ef m (k) = 2tcos(k y a) + 2t' [cos(2k x a) + cos(2k y a)] : (6) 

has its minima at the boundary of the reduced Brillouin 
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zone, k* = (±7r/2a, ±[tt — cos~ 1 (t/2t f )]/a) accounting 
for small electron pockets in MQO [3[, Fig.©. Placing 
the Fermi level near the bottom of the electron impurity 
band, \i = — 2.5t, yields the size of the electron Fermi 
pocket about 2% of the Brillouin zone as observed [3j- 
There is some anisotropy in the effective electron mass in 
the pocket, m x = H 2 /8t'a 2 and m y = h 2 /(%t' - 3t 2 /t')a 2 . 
Taking m x = 2m e [3|, and t' = 0.7t provides a reasonable 
estimate for t' ~ 32 meV and t ~ 46 meV. 

A momentum-distribution map of the impurity-band 
spectral function at the Fermi surface, A(k, 0) cx 
l/[(Ef m (k) - fi) 2 + 7 2 ] is shown in Fig.© with the in- 
verse quasiparticle life-time, 7 = 0.04t about the Din- 
gle temperature measured in MQO experiments Q. It 
hardly resembles the observed ARPES map [9], which is 
not surprising because the impurity-band ARPES matrix 
element, M(K) strongly depends on the photoelectron 
momentum K [17[. Calculating M(K) cx J <irexp(iK • 
r)ip(r) with the impurity-band wave function Eq.(pQ) one 
obtains M(K) cx fa , where /k is the Fourier component 
of the impurity-band Wannier orbital, f(~ft), introduced 
in Eq.(j3]). The extension of this orbital is comparable 
to the lattice constant, which explains the strong depen- 
dence of the matrix element on the in-plane photoelec- 
tron momentum Ky. 

Approximating f(~ft) as the 2D Coulomb bound state 
fim(p) yields 



Fermi surface, J(K, 0) cx M(K) 2 A(K,0), is 



M(K) cx 



1 



[l + (i^||a, m ) 2 ] 3 / 2 ' 
and the ARPES momentum-distribution map of the 



/(K ' 0)OC [l + (X || a im )2]3[( jB e m(K| 



)-M) 2 +7 



2T (8) 



A very close resemblance between theoretical Eq.(j8j) 
and experimental ARPES maps is shown in Fig.(j4j). The 
theory reproduces the locations, relative intensities and 
the half-moon shape of the ARPES spots in YBC06.5 
reconciling puzzling ARPES and MQO data. In fact, the 
small electron pockets observed in MQO experiments Q 
are seen in ARPES [9] , partially shadowed by the matrix 
element. The ARPES quasiparticle life-time turns out 
much shorter than in MQO, presumably due to an in- 
strumental broadening and a high level of impurity scat- 
tering off deposited potassium atoms on cleaved YBCO 
[9|. In the electron-doped cuprates, one expects an impu- 
rity band somewhat below the conduction band, if these 
compounds are doped charge-transfer insulators as the 
hole-doped ones. Then the outlined scenario with re- 
versed holes and electrons also accounts for the small hole 
pockets observed in MQO of the electron-doped cuprates 
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